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I. Introduction 


Controllability, observability, state feedtack, state estimation, and the 
notion of transfer matrix realization are all basically time domain concepts 
which have played an important role in the evolution of " modern control theory " 

; i.e. the analysis and design of dynamical systems via. state space techniques. . 
The primary purpose of this paper will be to present frequency domain analogies 
of these concepts in order to provide additional insight into the structure of 
linear systems, as well as offer alternative procedures for the analysis and design 
of linear multivariable systems. The reader should be aware at the outset that 
the results which will be presented here represent a natural extension of those 
presented in [I], Consequently, a clear understanding of the main ideas pre- 
sented in |lj would prove most helpful in reading this paper. 

The manipulation of polynomial matrices will also ploy a key role in the 
development employed. Therefore, several definitions and preliminary results 
regarding this class of matrices are given in section 2. In section 3, a basic 
result involving a frequency domain characterization of linear systems is proved. 
This result, which we will call a representation theorem for a rational transfer 
matrix employs polynomial matrices and is analogous to the time domain concept 
of a realization for a rational transfer matrix. The utility of this result, from 
the point of view of obtaining realizations, is then demonstrated in section 4. 

A practical application of transfer matrix representations is made in section 
5, which deals with the frequency domain analog of state feedback and estima- 
tion. Two examples are presented in order to clarify the various results presen- 
ted throughout the paper* and some concluding remarks are then made in section 6. 

2. Definitions and Preliminary Results 

A polynomial matrix , P(s), is any matrix whose elements ore polynomials 
( in this case in the Laplace variable s ). The degree of P(s) is defined as the 
degree of the polynomial or polynomials of highest degree comprising P(s). We 
will consider only finite degree polynomial matrices in this paper. 
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The following three elementary operations on P(s) may be defined: 

1) Interchange of rows ( columns ) i and | 

ii) Multiplication of row ( column ) i by a nonzero scalar 

iii) Replacement of row ( column ) i by itself plus any polynomial p(s) times 
any other row ( column ) j 

An elementary matrix , E(s), will be defined as any matrix which can be ob- 
tained from the identity matrix, I, by a finite number of elementary operations 
on I. Note that the determinant of any elementary matrix is therefore a non- 
zero scalar. Later we will show that ony polynomial matrix whose determinant 
is a nonzero scalar is an elementary matrix ( see Proposition 1 ). 

A polynomial matrix, P(s), will be called equivalent to another polynomial 
matrix, M(s), if and only if P(s) can be " reduced " to M(s) by a finite number 
of elementary operations on P(s); i.e. if and only if there are two elementary 
matrices, E|(s) and E^ts)/ such that E|(s)P(s)E 2 (s) = M(s). Equivalence will be 
represented by the symbol ^ ; i.e. P(s) — Mfc) . 

The rank of a polynomial matrix is defined as in the case of constant mat- 
rices; namely, the rank of P(s) is equal to the dimension of the largest minor of 
P(s) with nonzero determinant. A square polynomial matrix will be called non- 
singular if and only if its determinant is nonzero. The determinant of P(s) will 
be written as IP(s)l. The inverse of a nonsingular polynomial matrix will be 
written as P(s) or as P(s) + /IP(s)l/ where P(s) + represents the adjoint of P(s). 
Note that the only polynomial matrices, whose inverses are polynomial mat- 
rices, are elementary matrices. Furthermore, the adjoint of a polynomial mat- 
rix is a polynomial matrix. 

We now introduce the concepts of column proper ( and row proper ) poly- 
nomial matrices. In particular, let d^. denote the degree of the i-th column 
( an m x 1 matrix ) of the m x m matrix P(s). Then P(s) is of the form, 
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where the + ... denotes lower degree terms in s. Let be the constant m x m 

matrix ( p.. |, and diagl s^ c ' J the diagonal matrix with entries s^ c ‘. If P C (s) 

" ;l ^ ^ • 

is now defined as T diog( s Cl J, it can be verified by induction that lP(s)l = 
iris P + lower degree terms in s, where p = ]£d c |. ^ '* obvious that I P° (s)l = 
ys P , where “Y = II I. ^ 

P(s) will be called column proper if and only if the sealary / 0. Note 
that the " dual " concept of a row proper polynomial matrix can be defined in 
an analogous manner; i.e. P(s) is row proper if and only if its transpose, P^(s), 
is column proper. 

fVoposition 1 : Any ( m x m ) nonsingular polynomial motrix , P(s), can be re - 
duced to a column ( row) proper matrix , P(s), by a finite number of elementary 
column (row ) operations / 

Proof : If |P C (s)l / 0, P(s) would, by definition, be column proper and we would 
be done. On the other hand, suppose I P° (s) I - 0; i.e. 7= 0. This would imply 

that the m column vectors comprising P C (s) were Imearly dependent over the 

c c c 

monomials in s; i.e. if P.(s) represents the i-th column of P (s), then IP (s)l = 

0 would imply that 


i 

m 


£ P.(s)PfCs) = 0 
. 


( 1 ) 


for two or more nonzero monomials p.(s), i = 1 ,2, . . . . m. At least one of these 

monomials can be altered to unity by dividing (1 ) by a nonzero monomial of 

lowest degree, p, (s); i.e. ' 

m m 

£ Ip.W/p, MIP'M = £ p.M rf« - 0 (2 ) 

! • « 1 K I 1 

1=1 1 

Note that we need only establish this proposition for the case of column pro- 
per reduction - the proof for row proper reduction follows directly by duality. 
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where p^(s) = 1 . The replacement of column k of P°(s) by (2) Is analogous to 
postmultiplying P C (s) by the elementary matrix Ej(s), where Ej(s) is the identity 
matrix with an altered k-th column; i.e. 


E | ($) = 


10 ... 
0 1 ... 
. U 1 


P,(s) 0 . 
PjW 0 . 


pj) 0 

m 


0 

0 


(3) 


where at least one of the pf.(s), i / k, is nonzero. Pbstmultiplying P C (s) by 
E^(s) results in a zero k-th column of the product; i.e. if 


P cl (.) = P C (s)E,<»), 


C 1 • c 

then P (s) is equal to P (s) except for a zero k-th column. Now let 


P, (*) = PtoEjW 


(4) 


(5) 


Pj(s), thus defined, is equal to P(s) with an altered k-th column. In particular, 
the degree of the k-th column of P^(s) is strictly less than d^. Hie degree of 
the k-th column of P(s). Consequently, P^(s) is a new candidate for a column 
proper matrix; i.e. IP^(s)l =7^$^^ + lower degree terms in s, where pj is strictly 
less than p. It should now be obvious that if y^ = 0, this procedure can be re- 
peated as many times as necessary to produce ( for some { ) a nonzero y. . Then 
~ j ^ ^ j I 

P(s) = P(s)n E.(s) = P(s)E(s), where E(s) = II E. (s), thus establishing the proposition. 

Earlier in this section, we stated that any polynomial matrix whose determin- 
ant is a nonzero scalar is an elementary matrix. This fact can now be easily 
established using Proposition 1.- In particular, note that if IP(s)l were a nonzero 
scalar, then P(s) = P(s)E(s), where P(s) would, by the above arguments, be a 
constant matrix P. Therefore, P($) = Pf~ (s), and since any constant nonsingular 
matrix is an elementary matrix by definition, P(s), the product of two elemen- 
tary matrices, is itself an elementary matrix. 
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We now make use of the Smith canonical form tor polynomial matrices. 

I • i r. • . • i . t i 1 2 I . .I ii. i .1 


Employing a definition due to Rosenbrock, we will call two polynomial mat- 
rices, Pj(s) and Pjfc), of dimensions p x m and m x m respectively, relatively 
prime, if and only if the Smith canonical form of the composite ( p + m ) x m 


matrix 


m . P ml . 

r,(s)J IS [o J ; 


f and only if there exist two elementary matrices, 


E | (s) and £ 2 ( 5 ), of dimensions p + m and m respectively, such that: 


,(s) rpjsji E.(s) = rr 
KwJ 15 


A fundamental result pertaining to relatively prime polynomial matrices can now 
be established, namely: 


Proposition 2 : Let P (*) and P_(s) be relatively prime polynomial matrices of 


dimensions p x m and m x m respective! 


matrix M(s), there exist polynomial matrices, M.(s) and M_(s), of dimensions 


1 

m x p and m x m respectively, such that 

(s)P t (s) + M 2 (s)P 2 (s) = M(s) (7) 

Proof : We make direct use of the definition of relatively prime polynomial mat- 
rices; i.e. partition E^(s) of (6) into four submatrices, namely 


E|(s) = 


E„W E )2 (s) 
E 21 W E 22 (s) 


where E ^ ^ (s) is m x m, E^W is m x p, E^s) ispxm, and E 22 (s) is P x P • 
Equations (8) and (6) then imply that ( E^(s)P 2 (s) + Ej 2 (s)P^ (*) J = ^ m * 
Since E 2 (s) is an elementary matrix, we can postmultiply the above by the in- 
verse of E 2 (s) and premultiply by E 2 (s) to obtain: 

E 2 (s)E l 1 (s) P 2 (s) + E 2 (s)E 12 (s)P, (s) = I m (9) 

Premultiplying the above by M(s) establishes the desired result; i.e. equation 
(7), where 
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MWE 2 (.)£ I2 W«M,W (10) 

and 

M(i)E 2 (s)E ( ,(») = M 2 (») (II) 

3. Realizations and Representations 

The definitions and preliminary results on polynomial matrices established 
in the previous section will now be used to provide a rather natural analogy 
between time domain and frequency domain " factorizations " of linear systems. 
In particular, consider the following time domain ( state space ) description 
of a linear system: 

x = Ax + Bu ; y = Cx (12) 

where x is an n-vector, called the state, u is an m-vector, called the input, 
y is a p-vector, called the output, and A, B, and C are constant matrices of 
the appropriate dimensions. The triple { A, B,C | will be used as an alternate 
characterization of the system (12). If we take the Laplace transform of (12), 
and solve for y(s) in terms of u(s), assuming zero initial conditions on the state 
x, we readily obtain an expression for the transfer matrix ^ T(s), of the system 
in terms of the triple { A, B, C f ; i .e , 


y(s) = T(s)u(s), 

03) 

T(s) = C(sl - A)"'b 

(K) 


Note that T(s) is a p x m matrix of " proper " transfer functions; i.e. a proper 
transfer matrix . The term proper implies that each of the ( pm ) transfer func- 
tions comprising T(s) satisfies the condition that the degree of any numerator 
polynomial is strictly less than that of the corresponding denominator polynom- 
ial. If the degree of the numerator polynomial were allowed to be less than 
or equal to that of the denominator polynomial, the term realizable would be 

*We will assume throughout this paper that all possible pole-zero cancellations 
hove been made when speaking of transfer matrices. 
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used rather than proper. We will say more about real iza b le transfer matrices 
in section 5. 

Clearly, every system triple, j A, B, C J, yields a unique transfer matrix via. 
(14). However, the converse is not true; i.e. given any proper transfer matrix, 
T(s), one can readily characterize an entire class of triples which yield T(s) 
via. (14)1^ Any such triple will be colled a realization of T(s) of dimension • 
n, where n is the dimension of the A matrix. A realization of T(s) of lowest 

A . ^ A A . 

possible dimension n, denoted |A,B,C },wi II be called a minimal realization . 

The question of obtaining minimal realizations of proper transfer matrices 
has been the subject of numerous investigations[^'^^'^and it is well known 
that: 

i) Any Proper transfer matrix, T(s), does have o minimal realization, 

* r { A, B,C }, of dimension n. 

il) If { A, B,C | is a minimal realization of ^s). then the pair { A, B } is 
completely controllable, and the pair { A,c| is completely observable. 

The purpose of this section is to introduce a frequency domain concept 
which is analogous to the time domain concept of realizations of transfer mat- 
rices. In subsequent sections we will discuss and illustrate the utility of this 
concept. In particular, we define a controllable representation * of a p x m 
transfer matrix, T(s), as any pair,] (s), PjC*) },of polynomial matrices which 
satisfies the relationship: 

TW = P | (>)Pj'(.) (15) 

Clearly, Pj(s) and P^Cs) have dimensions p x m and m x m respectively. The 

order of a controllable representation will be defined as the degree of I P_ (s) I . 

,A a 4 

The pair,J Pj (s), P 2 (s) },will be called a minimal controllable representation of 

T(s) if and only if: 

i) The pair,{ Pj(s),P 2 (s) ), satisfies (15), and 

ii) ^ controllable representation of T(s) of order lower than that of 

{Pj(*)» ) ^ oes no * ex ' jt * 

In a dual fashion, we define an observoble representation of a p x m transfer 

*ln the next section, we will demonstrate how to obtain a controllable realization 
from a controllable representation. 
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matrix, T(s), as any pair, fQj(s),Q 2 (s) ] of polynomial matrices which satisfies 
the relationship: 

T« = Qj\»)Q 2 (.) (16) 

Clearly, Qj(s) and Q^s) have dimensions p x p and p x m respectively. The 
order of an observable representation will be defined as the degree of IQ.(s)l. 

A a ' 

The pair 'fQ^sKC^s) },will be called a minimal observable representation of 
T(i) if and only if: 

i) The pair,fQj(s),Q 2 (s) ), satisfies (16), and 

ii) on A observable representation of T (s) of order lower than that of 
f Qj(s),Q 2 (s)} does not exist. 

One additional definition is required prior to the establishment of a basic 

result of this paper. In particular, we recall that any proper p x m transfer 

aaa, a a 

matrix, T(s), has a minimal realization. (A, B,C J, of dimension n. If B is of 

A 

full rank m, and C of full rank p, the transfer matrix, T(s), will be called 
full rank proper . The ma|ority of physical systems with proper transfer mat- 
rices do satisfy this full rank condition, and most of our attention will focus 
on such systems. 

We will now establish a theorem basic to the structure of linear multivari- 
able systems and analogous to known results pertaining to time domain realiza- 
tions of transfer matrices. 

Theorem 1 : Any p x m full rank proper transfer matrix, T(s), has a m inima l 

A A A A 

controllable (observable ) representation , f Pj (s), Pj (») J ( {QjU^C^*)))* 
Furthermore, any such minimal representation satisfies the conditions : 

I) Jf (A,B, C ) is a minimal realization of T(s), then IPjWI ( IQj(s)l ) 
divides and is divided by Is! - Al . 

A A A A 

II) Pj (s) and P 2 ($) ( Q j (s) and Q^s) ) ore relatively prime polynomial 

matrices, 
t ' 

Proof : The proof of this theorem involves a combination of results. The first 

Vs in the case of Proposition 1, we need only establish this theorem for the 
case of controllable representations. Duality can then be employed. 


of these, the structure theorem, has been established elsewhere, and is repeated 
here for convenience as the first lemma; i.e. 

Lemma 1 : A ny p x m full rank proper transfer matrix , T(s), con be written os : 

T(.) = CSMI f ' (17) 

m 

A A.J 

where CS(s) and 6(») ore p x m and m x m polynomial matrices respectively . 
Furthermore, if (A,B,C ] is a minimal realization of T($), then IB !s(s)l = Isl-Al. 

A A m 

The proof of this lemma and the definitions of C, S(s),6(s), and B are given in 

in 

(I], and will not be repeated here. Furthermore, for convenience, we will let 

P,(s) = C$(s), and P (s) = B \(i); i.e. by (17) then, we can write T(s) = P.(s) 

*— 1 ^ m 1 

times (*)• 

The second result we will employ is, itself, a useful result, namely: 

Lemma 2 : fP^(s), P^(s) ) is on n-th order controllable representation of T(s), 

therr exists an n-th order controllable realization f A, B,C ) of T(s). 

A constructive proof of this result will be postponed until section 4, which deals 

with techniques for obtaining realizations fr«m representations. 

We note he-.e that Lemmas 1 and 2 can be employed to establish most of 

Theorem i . In particular. Lemma 1 establishes the fact that every full rank 

proper transfer matrix does have a controllable representation, namelyf P.(s), P_(s)), 

A a ' f 

and therefore a T! : ntmal controllable representation f P. (s), P_(s) J . We will, in 

fact, now establish that f P. (s),P.(s) ) is a minimal controllable representation. 

~ ^ a • A# A 

To do this, note that f P. (s), P_(s) ) is of order n, since IP-(s)l = I si - Al, where 
a ' * A l 

A determines a minimal realization of T(s) of order n. By Lemma 2, a controll- 
able representation of T(s) of order less than n can not exist, since this would 
imply the existence of a corresponding realization of T(s) of order less than n, 
which is impossible. Therefore, the pair f Pjfc)^^)) is a minimal controllable 
representation of T (s) . 

Lemma 3 : Consider the minimal controlloble representation, f P. (s), Pj(s) ) , of 
T(s) introduced in Lemma 1 . The Smith canonical form of 
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Proof : By Lemma 1 , I P^Cs)! - IsX ■ Al, where A is o member of a minimal real- 
ization of T(s). Using this fact, one need only repeat the arguments of Rosen- 
brock^ to establish Lemma 3. 

Lemmas 1, 2, and 3 establish the fact that given any p x m full rank proper 

transfer matrix, T(s), one can find a minimal controllable representation of T(s), 

namely { PjfsLPjU)} / which satisfies the two conditions of Theorem I . All 

that remains to be shown is that any minimal controllable representation of T(s) 

also satisfies these two conditions. We do this by establishing an equivalence 

relationship between minimal representations, namely: 

Lemma 4 : If {P, (s), P 0 (s)) is a controllable representation of the p x m full 
- j 2 ^ — — 

rank proper transfer matrix , T(s), and f P^ (s), Pjls)] , a minimal controllable 
representation of T(s), then there exists a nonsingular polynomial matrix , R(s), 
such that P^s) = Pj(s)R(s), and P^s) = P 2 (s)R(s); i.e. 

*_i 

(18) 


T« = PjMPj'w = P,(s)R(s)I P 2 (s)R(s) f' = P, (s)P^ ' (i) 


r 

Furthermore , IR(s)l divides I P^ts)! # and if f Pj (s), P 2 (s) ) is a minimal controll- 
able representation of T(s), R(s) is on elementary matrix* . 

V v.] 

Proof : Write T(s) as the product P. (s)P^ (s) as given by Lemma 1 . By Lemma 3, 


P| (s) and P 2 (s) are relatively prime. Hence, 


or 


TW=P 1 (s)P2 1 (s) = P 1 (s)P 2 , (s), 


P,(s) = 


P!(s)P 2 (s)P 2 (s) 

IP 2 (*)I 


(19) 

( 20 ) 


Since P^(s) is a polynomial matrix, IP 2 (s)l must divide the product Pj(s)P 2 (s)P 2 (s). 
Now, by Proposition 2, we can find a pair of polynomial matrices, fM^(s), 


M 2 (s) ) , such that: 


Mj (s)Pj (s) + M 2 (s)P 2 (s) =I m , 


(21) 


This final lemma is also important ir that it establishes the relationship between 
any controllable representation of a transfer matrix, and one of minimal order. 


or, 

M,WP,(»)p5wP 2 W * M 2 ( s )P 2 (s)P*( s )P 2 ( s ) - P 2 (s)P 2 (s) (22) 

Since P 0 (s)P 0 (s) - IP_(s)II , and since I P 0 (s)l divides the product P. (s)P_(s)P_(s), 
c l Z m i l it 

it follows that IP_(s)II can be factored from the left side of (22); thus IP 0 (s)l 
t m ^ z 

must also divide the right side of (22), namely "p (s) P^ (s) . Consequently, the 
quotient ( Pj (s) P 2 (s)/l P 2 ( s ) ^ ) ' 3 a polynomial matrix which we will call R(s), or 

Pj'WPjW =■ RW (23) 

Therefore, 

P 2 (s) = P 2 (s)R(s) (24) 

and from (20), 

P,(s)= P,(s)R(s) (25) 

It now follows from (24) that IPjWl = IPjWMRMl, or that IR(s)l does divide 
IP^ts)! . Furthermore, note that iff PjfebPjfs)) were a minimal controllable 
-representation of T(s), then I P^Cs)! would be a polynomial of degree n, as would 
I P 2 (s) I • Since we have just established the fact that ^(s)! = I Pj (s) IX|R (s) I , 
this implies that IR(s)l must be a polynomial of degree zero; i.e. a scalar. 
Consequently, if f P^ (s), Pjt 5 ) ) is a minimal controllable representation of 
T(i), R(s) must be an elementary matrix. Lemma 4 is thus established. 

• Therefore, since f P.(s),P 0 (s) ) satisfies the two conditions of Theorem 1, any 

A A ^ 2 

{ = f P| (s)E(s) # P 2 (s)E(s) ) also satisfies these conditions by direct 

arguments. Theorem I is thus established . 

To summarize thus far, we have introduced a new concept pertaining to 
the transfer matrix characterization of linear systems. This frequency domain 
concept, called a representation, is analogous to the time domain concept 
of realization. More specifically, we have shown that any full rank proper 
transfer matrix, T(s), can be factored in three different " minimal " ways; i.e. 

A A A A A.l A _] A 

T(s) - C(sl - A) B = P, (s)P 2 (s) = Q , (s)Q 2 (s) (26) 

All of these " factorizations " share certain common properties. In particular. 
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A A A A 

IsI-AI, I (s) I # and IQ ^ (s) I are all n-th degree polynomials differing by, at 
most, a scalar multiplier. The word " minimal " associated with these three 

A A - ] A 

factorizations implies ( as far as the quantity C(sl - A) B is concerned) that 

A .A A A 

the pair f A, B ) is controllable, and the pair f A,C J is observable. An onal- 
ogous interpretation regarding the pairs f Pj^P^s)] and f Q ^ (s), Q^Cs) J will 
be the subject of the next section, which deals with the derivation of realiza- 
tions from representations. 

4. The Derivation of Realizations from Representations 

The results which have been presented thus far are primarily of academic 
interest; i.e. if one could obtain a controllable or observable representation 
of some transfer matrix, T(s), it is not at all clear how such a representation 
could then be used for the purpose of analysis or design. The remainder of 
this paper will focus on such practical considerations. Before presenting any 
direct utilization of representations however, it will be of interest to consider 
the relationship between realizations and representations. 

For reasons which will become more obvious as we progress, we will first 
consider the question of converting from a representation* to a realization. 

In particular, we will demonstrate constructively how one can always obtain 
an n-th order realization from an n-th order representation, thus establishing 
Lemma 2 of the previous section . 

Consider the n-th order representation, f } , of the full rank 

proper transfer matrix T(s); i.e. IP 2 (s)l is a polynomial of degree n. By Prop- 
osition 1 ( section 2 ), P^^) can be reduced to column proper form via. some 
elementary matrix E*(s). Call the resulting column proper matrix P^s); i.e. 

PJ(s)= P 2 (s)E*(s) -- (27) 

Let Pf(s) = Pj(s)E*(s), and note that 

TW-PfWFjf'w (28) 

* Again, we will deal primarily with controllable representations— by duality, 

% 

analogous results hold for the case of observable representations. 
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f P|(s)/P*( s ) 3 * 5 also a control lablr n:presr ntation of T(s) of order n. 

Define G* as the constant m x m matrix consisting of the coefficients of 

the highest degree s term or terms in each column of P*(s) ( G* is analogous 

to the matrix Fdefined in section 2). G* is nonsingular since IPJ C (s)l is a 

-1 * 

nonzero monomial. If P£(s) is premultiplied by G* , the resulting polynomial 
matrix, which we will call D(s), is characterized by the fact that the highest 
degree s term in each ( i-th ) column of D(s) appears in the corresponding ( i- 
th ) row with unity coefficient; i.e. 

D(.) = G*‘ I PJ( S ) (29) 

where D(s) can also be written os: 

D(s) = D c (s) - A $(s), (30) 

m 

where D C (s) is a diagonal matrix with entries s^' , i = 1, 2, ... m. The trj 
have a special interpretation which is more thoroughly discussed in [1J. A 

m 

is an m x n constant matrix whose elements are directly determined by (29) 
and (30), and S(s) is an n x m matrix of monic monomials in s. In particular, 



We can write T(s) as the product: 


TM - PJMI G*G*' , P‘M )■' =• P*MD‘ 1 MG*' 1 . (32) 

This expression for T(s) admits a direct time domain realization ( A, B,C ) by 
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virtue ot the results given in jl |. In porticulor, P*(s) is analogous to CS(s) , 

[1| _] A I] I 1 

D(s) to6(s), and G* to B 

m 

More specifically, an n-^h order realization, f A,B,C ] , of T(s) can be 

obtained as follows; An appropriate p x n matrix C can be found directly by 

inspection, since PT(s) can be written as CS(s)^j where S(s) is given by (31). 
k 

Now define d, = as in (1J. Note that d n. Replace the ( m ) d, -th 


1 


m 


rows of an n dimensional companion matrix by the ordered ( m ) rows of A , 


m 


where A is given directly by (30) in much the same way as C was obtained. 

m . . [1J 

The resulting matrix is an A corresponding to the above choice for C . Sim- 
ilarly, replace the ( m ) d. -th rows of an n x m null matrix by the ordered 

k . -1 

( m ) rows of the m x m nonsingular matrix G* . The resulting matrix is a B 
corresponding to the A ond C already selected. The triple f A,B,C ] , thus 
obtained, is a controllable realization of T(s) obtained from the controllable 

representation, { P,(s),P 0 (s) ) , of T(s). Herein lies the motivation for associat- 

1 i m -1 
ing the word " controllable " with a representation which factors as PjfcJPj (*)• 

By duality, one could obtain an observable realization of T(s), given an ob- 
servable representation, f Q j (s^C^U) ) , of T(s). A simple algorithm to use 
is the one outlined above, only applied to the transpose of T(s), and later 
transposed again; i.e. T^(s) = would yield an observable triple 

( realization ), { aJ cj B^ ) . 

The algorithm which we have just outlined for obtaining realizations from 
representations is just that, an algorithm. The reader should really be com- 
pletely familiar with the results given in [1] in order to interpret it as a con- 
structive proof of Lemma 4. 

So far, we have discussed certain implications associated with certain 
representations of T(s). However, we have not discussed a more fundamental 
question, namely how to obtain representations; i.e. factorizations of T(s) 
as either the product P^sJPj^s) or Q ^ ^ (sJQjW • Also, we have not presented 
any methods for determining whether or not a particular representation is min- 
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imal. Some discussion of these questions appears to be in order. 

First, we remark that in general, the problem of finding a representation 
of T(s) is analogous to that of finding a realization of T(s). In fact, the al- 
gorithm just outlined allows one to rather easily obtain an n-th order realiza- 
tion from an n-th order representation of T(s). The converse is also straight- 
forward and is one of the results presented in [1], We might also remark that 
it is relatively easy to obtain either a realization or a representation of a given 
transfer matrix provided minimality is not important. Unfortunately, in most 
applications, minimality is a key concept, since it is tied to controllability 
and observability of the system under consideration. In the case of realizations, 
" reducing " any n-th order realization to an n-th ( minimal ) order one in- 
volves some rather intricate computations!^ An analogous condition holds in 
the case of representations as one might expect; i.e. as previously shown - 
see Lemma 4 - if a controllable representation f P^sJ^^ts) ) is not minimal, 
there exists a nonsingular polynomial matrix, R(s), which can be postmuiti plied 
( factored ) out of both P^(s) and P 2 (s) to produce a minimal controllable rep- 
resentation. 

The development of algorithms for determining when two polynomial mat- 
rices are relatively prime, finding postfactors such as R(s), and factoring poly- 
nomial matrices has and will continue to be the subject of further investiga- 
tions. Some interesting results have already been obtained. 

5. Feedback Compensation 


In this section, we will present an application of the methodology develop- 
ed in the previous sections. In particular, we will outline a ( frequency domain ) 
feedback compensation technique for linear systems in order to achieve some 
desired closed loop transfer matrix. An analogy between this technique and 
( time domain ) state estimation and feedback will then be made and demon- 
strated by example. 
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mmmm - * 


To begin, consider an/ controllable rupiesentation, f P ^ (s) ^ P2 (s) } , of 
T(s), where (s) is ossumai ( for convenience ) to be column proper and equal 
to G*D(s) ( see eqs. (29) and (30) ). Suppose we desire the " closed loop " 
transfer matrix Tp(s), where 

T f (s) = P,(«)P‘p(s), (33) 

and 

P 2F (») = P 2 (s) - FS(s), (34) 

S(s) is given b/ (31), and F is an m x n constant matrix. We will call the 

pair f P.(s),P-_(s) ) , a closed loop controllable representation of the transfer 

’. t ZP 
matrix . 

Let us now consider a design procedure which will produce this desired 
closed loop transfer matrix under the appropriate conditions. In particular, 
consider the following block diagram, which represents feedback compensation 
for the system whose open loop transfer matrix T(s) = P ^ (s) Pj ^ (s); i.e. 

y(0 


Figure 1 

Frequency Domain Estimation 

We define a frequency domain estimator of FS(s) - PjW ~ PjpW* os o 
triple. f Qjts^G^s^Q^U) ) / of polynomial matrices of dimensions m x p, 
m x m, and m x m respectively, which satisfies the following three conditions: 

*The motivation for this terminology will become more apparent when we later 
relate F to closed loop " state " feedback in the time domain. 
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(i) Q 3 (s) is nonsingular, and IQ^s)! is Hurwitz; i.<*. the roots of IQ^s)! 
lie in the half-plane Ke s< 0. 

(ii) Q 3 \s)( (s) # Q ^ (s) J is a realizable ( section 3 ) m x ( m*p ) tionsfer 

matrix; relating the Laplace transform of the output, v(s), of the es- 

u(sf 


timator, to the Laplace transform of the input 
(iii) The following relationship holds: 


y(s) 


Q 3 (s)FS(s)= Q , (s) P | (*) + Q 2 (s)P 2 (s) 


(35) 


The following important fact can now be established: 

Theorem 2 : The above three conditions which define a frequency domain es - 

timator are sufficient to insure that the closed loop transfer matrix, T (s), de- 

- 1 ^ 

picted in Figure 1 ; i.e. y(s) = Tp(s)w(s), is equal to Pj(s)P 2 p(s) os defined by 
eqs. (33) ggd (34). 


Proof : Consider the equations implicit in Figure 1; i.e. 


Also, 


u(s)= w(s) + Q 3 '(s)Q 2 (s)u(s) + Q 3 \s)Qj(s)y(s) 


-1 


y(s)= P 1 (s)P 2 '(s)u(s) 


Substituting (37) for y(s) in (36) and collecting terms, we obtain: 

II- Q 3 ' (*) l Q 2 (0 ■ QjMPjMP^W 1 1 u (*) = w (*)* 


or 


Q*' (0 [ Q 3 (s)P 2 W - Q 2 (s)P 2 ( s ) - Q ] (OP, (0 I Pj' (.)«(») = w(s) 
If (35) is now used in (39), we obtain: 


I P 2 (s) - FS(s) J P^frMs) = w(s). 


or employing (33), 




Since P 2 p(s) '* non *' n gular. 


ID 


-1 


u(») = P 2 (s)P 2 ;(sMs), 


(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 
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( 43 ) 


and using (37), we obtain the desired result; i.e. 

y(s) = P 1 (s)P 2 p(s)w(s) 

which establishes the theorem. 

To summarize, we have now shown how one can design a feedback compen- 
sation system (a frequency domain estimator ) which produces a desired closed 
loop transfer matrix. Note that no specific reference to the " state " of the 
system has been made, nor has it been necessary. However, it is perhaps app- 
ropriate to now employ this concept in order to relate our results to known time 
domain results. In particular, there is a rather natural analog between the com- 
pensation scheme just presented and the time domain concepts of state estimation 
and feedback, as we will now show. 

To begin, let us employ the algorithm given in section 4 to obtain a contro- 
llable realization f A,B,C ) , of T (s) corresponding to f Pj(s),P 2 (s) J ; i.e. 

T(s) = PjMPj'w = C(sl - A) -, B (44) 

Recalling the interpretation of the realization f A,B, C ) in the time domain 
via. the differential equation (12), we will define state feedback ( in the time 
domain ) as the control law, 

u = Fx + w, (45) 

where F is the m x n constant matrix given by (34), and w is an m-vector rep- 
resenting the external input to the system as indicated in Figure 2; i.e. 



Figure 2 

Closed Loop State Space System 
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Substituting (45) for u in (12) results in the closed loop system, 

x = (A ■» BF)x + Bw ; y ~ Cx (46) 

as represented by Figure 2. 

The closed loop transfer matrix associated with this system is readily deter- 
mined; i.e. 

y(0 = T p (s)u(s), (47) 

where 

T p (s) = C(s I - A - BF)" 1 B (48) 

Equation (48) clearly contains a rather simple alteration of the original expres- 
sion for the open loop transfer matrix as given by (14), namely the replacement 
of A by A + BF. The analogous interpretation of " state feedback " from the 
point of view of the controllable representation, f , of T(s) has 

already been made, namely the replacement of P^(s) by Pjpfs) 05 9' ven by (34). 
In short, from our analysis thus far and some rather elementary arguments given 
in [1], it follows that: 

T f « = T f <0 = C(il - A - BF)' 'b = PjWPjJw (49) 

State feedback in the time domain can therefore be directly related to a closed 

loop controllable representation of the transfer matrix. 

State estimation in the time domain can similarly be related to frequency 

domain estimation. Due to spatial limitations however, we will not present 

a time domain characterization of state estimation here since we would only 

(81 [91 

be repeating well known results. We should emphasize however, that in 
order to achieve a state feedback design ( as depicted in Figure 2 ) when only 
the output y of the system is measurable, it is necessary to employ a technique 
such as state estimation or frequency domain estimation. 

Note that we have thus far avoided any reference to the term " minimal M 
when discussing either realizations or representations of the systems under con- 
sideration in this section. It is, nevertheless, an important concept for the 
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following reason. In the case of time domain state estimation, it is well known 

that in order to construct a total state estimator, s/stem observability is requir- 

( 8 ] 

ed; i.e. the controllable realization f A,B,C )we have been discussing, must 
also be observable, and hence minimal Similarly, from the point of view 
of frequency domain estimation, a solution to (35) is guaranteed, regardless 
of the choice of Q^(s), provided the representation f P^s^PjU) ) is minimal; 
i.e. Theorem 1 then guarantees that Pj(s) and P^(s) are relatively prime and 
hence that Proposition 2 holds. 

A number of unresolved questions associated with this frequency domain 
approach to design remain which will require additional investigations. For 
example, nothing has been said regarding the increase in system order required 
to implement the frequency domain estimator. It might be desirable to keep 
the degree of IQ^(s)l ( the order of the estimator ) as low as possible. The prob- 
lem of finding a triple ,fQ|(s),Q 2 (s), Q^s) ), which produces the desired closed 

loop transfer matrix is not a trivial task in the case of multivariable systems ( 

(81(91 

the existence of such a triple depends on known time domain results 1 ) . 
Unfortunately,. we cannot answer all questions related to the implementation 
of frequency domain estimator feedback at this time. Under certain conditions, 
however, the construction and feedback utilization of frequency domain estim- 
ators is relatively simple and straightforward, and we will focus our attention 
on two such cases in this paper. In particular, 

Cose 1 ) Consider the case when: 

a) Pj(s) is square and nonsingular 

b) IP (s)l is a Hurwitz polynomial 

-1 . 

c) FS(s)Pj (s) is a realizable transfer matrix 


If in (36), (^(s) is set equal to zero and 


FS(s)P'\s) * Q^WQjts), 


(50) 


then substituting FS(s)P^(s) for Q^(s)Qj(s) in (36) and employing (37) directly 
yields the relationship: 
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I'W *(•>) • » S(s)p^ 


II - FS(i)P" (s)|u(s) w(s), (52) 

from which equations (40) through (43) follow directly, thereby producing the 
desired closed loop transfer matrix T^(s) P^P^s). 


Example of Cose 1 : 
Suppose that 


[ s 3 - 3s 2 - I Os * 3, s 3 +3s~| 
-3s 3 - s 2 ♦ 3s - 17, 2s 2 ♦ 6j 


4 3 2 

s ♦ s - 5s * 3s 


The reoder can verify that a minimal controllable representation, f PjfshPjM ) r 
of T(s) is the pair, 


P,(s) - 


~s *■ 3 i r 

-3s - 7, 2J ' 


where 


IPj(s)l - 3s 2 f 9s ♦ 6 -• 3(s ♦ l)(s ♦ 2), 


( note that conditions a and b defining Case 1 are thus satisfied ) 
and - 


where 


Remark: In this example, P,(s) and P 0 (s) con be obtained using the algorithm 

17| 1 

of Menahem to factor the numerator of T (s) . 

Note that this open loop system is clearly unstable. Suppose we now wished an 

asymptotically stable N decoupled " closed loop system. Employing the results 

A 

given in ( 1 J and (10], we conclude that in this case, a Pjp(s) and G exist such 
that: 


A 

“> 2 o,. 0 1 


P 2 (») ^ 

_4s 4 5 s 2 - 2s + lj ' 

(E4) 

iPjMi * 

s* ♦ s 3 -5s 2 ♦ 3s = s(s ♦ 3)(s - 1) 

(E5) 

A 

A 
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M 


A A_i 

P l ls)P 2 F (s)G 


s * 4 ' 

0 , 


after oil possible pole-zero cancellations have been made. We directly deter- 

. #fc . iiojiii] 

mine that : r- — » 


0 1 

1 -I 


Since P^(s) and G are now known, we can solve (E6) for P 2 p(s) = G 


times Pj(s) # and for this example. 


s + 4, 0 

~ s + 5 


^ s 2 + 22 / 3 * + 35/3, 

F> 2F (s) ' -l/3s + 1/3 , s 2 + 


-2/3s - 10/3 
14/3s + 10/3 


FS(s) is now given by (34); i.e. 

FS(s) = P 2 (s) - P 2p (s) = 


-22/3s - 26/3, 2/3s + 10/3 
_13/3s + 14/3, -20/3s - 7/3 


Solving for FS(s)P (s) using (E2) and (E9), we obtain: 

T~ 2 2 

I A * . J ft ^ . 1 


FS(s)p7'(s) 


2s 2 + 6 , 8s+14s+10~ 

= L“20s 2 -45s -7, -11s 2 -27s -7_ ■ 


(E10) 


3(s + l)(s + 2) 


Clearly, FS(s)Pj (s) is a realizable transfer matrix( condition c of Case 1 is 
thus satisfied ), which , when substituted into the feedback path as indicated 
below in Figure El, yields the desired closed loop transfer matrix as given by 
(E6). 



Figure El 

Cose 1 Closed Loop System 
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The overall system represented by Figure Cl is 6-*b order with four pole-zero 

cancellations at s r -I, -1,-2, and -2. An algorithm for constructing a time 

A_ | 

domain realization of F5(s)Pj (s) can easily be determined using the results 


presented here ( section 4 ) and in inference (I 

z = Hz + Jy 
v = Kz * Ly 

where 


In particular, if 


(EH) 


i-7/3 2 /ji To - 1 / 3 ] f 9 -*] 1 . rv> 3/3 

H ~ -2/3 -2/'3 ' J " I 1/3 ' K ~ -1 5 ' L= -20/3 -11/3 ' 


then: 


v-1 




K(sl - H) J U = rS(s)Fj (s) 


(£ 12 ) 


as given by (blO). 

Case 2) Consider the case when m = p = 1; i.e. the case of scalar ( single in- 
put, single output ) systems. In this case, our results can be tied to th» classical 
resultant and eliminant matrix of Sylvester. In particular, we will write the open 
loop transfer function T(s) as: 


T(0 = 


P,(s) 

P 2 (s) 


p io + p n s + •••* 


+ P l,n-1 $ 


n-1 


(53) 


p 20 +p 21 S+ ”• 


** +p 2n S 


The eliminant associated with this transfer function (or the pair of polynomials 
f Pj( s ),P 2 ( $ ) ) ) is the ( 2n - 1 ) x ( 2n - 1 ) matrix E defined below; i.e. 


(n- 1/ tows 

E = |0 ... 0 p in p,, ... p, ,| v (54) 

n rows 


0 ... 0 

p 20 P 21 •“ 

" P2n ll 

0 • • • P 2 o 

P 21 *•* 

•; P 2n ° 

• • 


• i 1 

m I 

p 20 P 21 • • 

P 2n 

0 ... 0 / 

0 ... 

0 p 10 P 11 

P 1,n-I \ 

0 ... 

. . • 

p io. p n 

P l,n-I ° 

• 9 1 

• • 
• 

• • 

• f 

p 10 P 11 •• 

. p, , 0 
r l,n-l 

0 1 


*A discussion of Sylvester's resultant and eliminant matrix con be found in a 
number of earlier texts dealing with linear equations; e.g. (12). 


Note that the entries of E, the eliminanf, are simply the coefficients of the 
polynomials Pj(s) and p 2 (s). Sylvester first defined the eliminant matrix and 
its determinant , the resultant R, and established the fact that the resultant 
R = IEI is nonzero if and only if (s) and P^s) are relatively prime . The sig- 
nificance of thi< fact, from the point of view of frequency domain estimation, 
is most important as we will now show. 

In the scalar case, a frequency domain estimator is o triple, f Q ^(i) , Q i 
Q^(s) ) , of polynomials which satisfies the three conditions discussed earlier 
in this section. In particular, Q (s) is a polynomial whose roots are the "poles" 
of the estimator. It is well known that in the scalar case, complete observa- 
bility of a system is sufficient to insure that a ( time domain ) estimator with 
arbitrary poles of order n-1 can be constructed . 1 Furthermore, this estima- 

tor is characterized by the fact that its state and the system output y "exponen- 
tially approach " the state, x, of the open loop system with time. An analogous 
result holds in the frequency domain. In particular, let Q-j(s) be any arbitrary 
Hurwitz polynomial of degree n-1. Recall that FS(s) , the difference between 
P^ts) and Pjp (s) is a polynomial of degree 4 n-1 . Therefore, the product Q^s) 
times FS(s) , as given by (35), is a polynomial of degree 4 2n-2. This product 
will be written as: 

2n-2 

Q 3 (s)FS(s) = m Q + m j s + ... ,, + m 2n-2 S (55) 

The right side of (35) must also be a polynomial of degree 4 2n-2, and since 
Pj(s) and P^(s) are polynomials of degree4n-l and n respectively, it follows 
that Qj(s) and Q^is) must be of degree no greater than n-1 and n-2 respect- 


ively; i.e. 

QjW = ^10 + q u s + *•* 

n-1 

i, + q l,n-l* ' 

(56) 

and 

Q 2 <‘> = q 2 o + q 2i s f 

n-2 

** + q 2,n-2 S 

(57) 


Note that this also insures that (s) ( (^(shQjfc) J is a realizable transfer 
function. The reader can verify that the right side of (35), namely Qj(s)Pj(s) 
+ is given by: 


U 

0 

s 

0 

0 

I 
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I 

B 


lq 2,n-2' q 2,n-3' “* q 20' q I ,n- 1 ' * * * q 10 jE 


2n-2 


I *6 • 


QjWPjW+QjWPjW * 


(58) 


Q^(s)FS(s) # as given by (55), can now be written more succinctly as m*s*, 

where m* = ( m^, m^, m 2n-2^‘ ^*‘ n 9 ^‘ s no * a *' on ' (35) reduces 

to; 

m*s* = q*Es # (59) 

Clearly, a solution, q*, of (59) exists if E is nonsingular, namely 

q* = m*E _1 


(60) 
t 


Equation (60) can be employed whenever Pj(s) and (s) are relatively prime. 

Remark: Equation (59) has been extended to the multivariable case ( for terrain 

‘ [I3j 

classes of systems ) by the author. 

The utilization of equations (59) and (60) for feedback design will now be ill- 
ustrated ; i.e. 

Example of Case 2 : 

2 

Consider the open loop system whose transfer function T(s) = (s+3)/(s +s-2); i.e. 

2 

(s + s -2 )y(s) = (s + 3) u(s). Clearly, P^s) = P lQ + p^s, where p 1Q = 3, and 

p ] j=l . Also, P 2 (s) = p 2Q + P 2] s + p 22 $ 2 , where p 2Q = -2, p 2] = 1 , and p 22 = 

1 . Furthermore, in this case, Pj(s) and P 2 (s) are relatively prime, as the reader 

can readily verify; i.e. P 2 (s) = (s - 1 )(s + 2). The open loop system is thus 

unstable. Suppose we wish to design a feedback control system which would 

yield closed loop poles at s = -1 and s = -4. This would imply that P (s) = 

2 

4 + 5s + s ; i.e. FS(s) = -6 -4s, and if we arbitrarily choose Qg( s ) = s + 6, 
then Q 2 (s)FS(s) = -36 -30s -4s 2 , or m* = ( -36, -30, -4 J ond s* = i in 
accordance with (55), (58), and (59). By (54) and (59), 

*This condition is equivalent to complete controllability and observability of 


the open loop system ( in the time domain ). 


(14J 
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1 


q 2 Q* q ir q io 

q* 


0 3 1 


3 1 ol s 2 


and by (60), 

q* = m*E‘' * | 9/2, -17/2, -9 ] ( , 

C q* ( Qj(0 and Qjfa) ) ^ as ^ een determined, the closed loop system con 
be constructed as indicated by Figure 1. For this example, one can readily 
verify that the transfer function of the closed loop system depicted below is 
equal to (s + 3)/(s + 1 )(s + 4) as desired. The system shown is actually a third 
order system with a pole-zero cancellation at the pole of the estimator; i.e. 




Figure E2 

Closed Loop Scalar S/stem 


One final point is worth noting before we conclude, namely if in (36), we 
let Qj(s) = 0 and Q^VjG^s) = FS(s) ^ (s), we would then have 


u(s) = FS(s)P 2 (s)u(s) + w(s). 


( I- FS(s)P’ (s)lu(s) -- w(s) 


from which equations (40) through (43) follow directly; i.e. the control law rep- 
resented by (61) yields the dnvred " closed loop " transfer matrix T_(s) ~ P.(s) 

-1 r 
times P 2 p(s) without 0 7 o •*» * feedback, since Qj(s) ~ 0 ( see Figure 1 ). This 

input feedback " control low represented by (61) can also be represented by 






k 


n 


1 
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li 


J 





the relation: 


( 63 ; 


'j(') - P 2 (*) P 2F^ W ^ ; 

i.e. from (62), I P^s) - FS(s) |P 2 '(s)u(s) - w(s), and since P^(s) - FS(s) = 
P 2p (s), (63) follows directly. No^e that (63) simply represents feedforward 
compensation. Whot this simply means is that any " state 11 feedback control 
low can be realized by feedforward compensation . Note however, that the 
control low represented by (61) (or (63) ) implies: 

i) pol e-zero cancellations at the zeros of I P^ (*) 1 # and 

ii) an increase in system order equal to n, the order of the given system. 
Neither of these two conditions is desirable in most cases for rather obvious 
reasons; i.e. if I P^ (s)l were not o Hurwitz polynomial, the feedforward com- 
pensated system would be unstable. Also, if one employs o feedback design, 

[8] (9) 

the increase in system order can always be kepti n-m. J There are other 
reasons for employing feedback rather than feedforward compensation in most 
cases. A primary reason is sensitivity reduction, a subject which will not b: 
discussed in this pcper, but one which is quite important, especially when con- 
sidered from the point of view of linear optimal control 
6. Concluding Remarks 


We have introduced a new concept and demonstrated its utility in the analysis 
and design of linear systems. This ( frequency domain ) concept, namely the no- 
tion of representation for proper transfer matrices, was shown to be analogous to 
the ( time domain ) concept of realization. An algorithm for obtaining realiza- 
tions from representations was presented. Frequency domain analogies of other 
time domain concepts were also given, notably a frequency domain character- 
ization of ( simultaneous ) state estimation and feedback. For the most part, the 
results presented relied heavily on some rather elementary definitions and results 
dealing with polynomial matrices and their manipulation. Areas for future work, 
tying a closer bond between polynomial matrices and linear system theory, were 

pointed out at various times. Additional results employing the techniques dev- 

, (131(15] 

eloped here will also appear soon. 
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